Self-Regulation of Solar Coronal Heating Process via Collisionless Reconnection 
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I propose a new paradigm for solar coronal heating viewed as a self-regulating process keeping 
the plasma marginally collisionless. The mechanism is based on the coupling between two effects. 
First, coronal density controls the plasma collisionality and hence the transition between the slow 
collisional Sweet-Parker and the fast collisionless reconnection regimes. In turn, coronal energy 
release leads to chromospheric evaporation, increasing the density and thus inhibiting subsequent 
reconnection of the newly-reconnected loops. As a result, statistically, the density fluctuates around 
some critical level, comparable to that observed in the corona. In the long run, coronal heating can 
be represented by repeating cycles of fast reconnection events (nano- flares), evaporation episodes, 
and long periods of slow magnetic stress build-up and radiative cooling of the coronal plasma. 

PACS numbers: 96.60.lv, 52.35.Vd, 94.30.cp, 96.60.P- 



This paper is devoted to the problem of solar coronal 
heating (see Ref. [l| for a recent review), viewed in the 
context of Parker's nano-flare model 0]. Since the main 
heating process in this model is magnetic reconnection, 
I will first summarize the recent progress in reconnec- 
tion research achieved in the past 20 years. Even though 
a complete picture of reconnection is still not available, 
there is now consensus about some of its most important 
aspects. My main goal is to apply this new knowledge to 
the coronal heating problem. 

First, I want to emphasize the importance of a realiza- 
tion by Petschek 0] that the main bottleneck in the clas- 
sical Sweet-Parker 0, [B| reconnection model is the need 
to have a reconnection layer that is both thin enough for 
the resistivity to be important and thick enough for the 
plasma to be able to flow out. Furthermore, Petschek [3[ 
proposed that this difficulty can be mitigated if the recon- 
nection region has a certain sub-structure — the Petschek 
configuration, with four shocks attached to a central dif- 
fusion region. This results in an additional geometric fac- 
tor leading to faster reconnection. This idea is especially 
important for astrophysical systems, including the solar 
corona, irrespective of the actual microphysics inside the 
layer. Indeed, the system size L is usually much larger 
than any microscopic physical scale S, e.g., the ion gyro- 
radius pi, the ion collisionless skin-depth di = c/u>pi, or 
the Sweet-Parker layer thickness <5sp = ^Ltj/Va- There- 
fore, a simple Sweet-Parker-like analysis would give a re- 
connection rate v vcc /Va scaling as 5/L <C 1, and hence 
would not be rapid enough to be of any practical interest. 
Thus, we come to a conclusion that Petschek's mecha- 
nism (or a variation thereof) is necessary for sufficiently 
fast large-scale reconnection. 

Recently, however, several numerical and analytical 
studies (e.g., @, 0, 0, ©, El E3, E3, E3 El El, El) and 
laboratory experiments [1 71 ] have shown that, in resis- 
tive MHD with a uniform (and, by inference, Spitzer) re- 



sistivity, Petschek's mechanism does not work; the slow 
Sweet-Parker scaling applies instead. In other words, In 
the collisional regime, when classical resistive MHD ap- 
plies, one does not get Petschek's fast reconnection. 

It is then natural to ask now is whether fast reconnec- 
tion is possible in a collisionless plasma where resistive 
MHD is not valid. The answer now appears to be "yes" . 
First, in space and solar physics fast collisionless recon- 
nection events have been observed for a long time. More 
recently, it has also been seen in laboratory experiments 
[171 Il8|. In addition, several theoretical and numerical 
studies have recently indicated that fast Petschek-like 
reconnection does indeed take place in the collisionless 
regime. Moreover, there may even be two physically- 
distinct mechanisms for fast collisionless reconnection: 
(1) Hall effect (e.g., 0, El, El ED, 0, Hi); 
and (2) s pat ially-localized anomalous resistivity (e.g., 
0, !, [ll, Hi, EE El, El, [H, (13). At present, it is 



still not clear which of them operates under what con- 
ditions and how they interact with each other. How- 
ever, both of these mechanisms seem to work and both 
seem to involve an enhancement due to a Petschek-like 
configuration. Thus, I believe it is safe to say that a 
Petschek- enhanced fast reconnection does indeed happen 
in the collisionless regime. 

To sum up, there are two regimes of magnetic recon- 
nection: the slow Sweet-Parker reconnection in resistive- 
MHD with classical collisional resistivity, and the fast 
Petschek-like collisionless reconnection. 

Now, how can one quantify the transition between 
these two regimes? First, consider the case with a rel- 
atively weak (or zero) guide field, B z < B , where B 
is the reconnecting field component. Then, the condition 
for fast collisionless reconnection can be formulated (e.g., 
d, El, El 0, 0, [H ) roughly as 



S S p < di . 



(1) 
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(Since the discussion in this paper is very approximate, I 
will consistently ignore all numerical factors of order 1.) 

Expressing resistivity that enters via Ssp in terms of 
the Coulomb-collision electron mean- free path A eim f p , one 
gets 0: 
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where (3 is the ratio of the plasma thermal pressure 
(2n e T e ) at the center of the layer to the reconnecting 
magnetic field pressure (B 2 /8ir) outside of the layer. The 
condition of force balance across the layer (in the ab- 
sence of a strong guide field) requires /3 ~ 1, where we 
neglected the contribution due to the upstream gas pres- 
sure. Then, the above fast collisionless reconnection con- 
dition becomes 



L < L c = \J nii/m e A 
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40 A 
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Note that, by construction, the mean- free path that en- 
ters here is due to classical Coulomb collisions. It is given 
by A e ,mfp — 7 • 10 7 cmnj" 1 T 2 , where we took log A ~ 20 
and where n 10 and T 7 are the central layer density n e 
and temperature T e in units of 10 10 cm -3 and 10 7 K, 
respectively. Substituting this into Eq. we get 



L < L c {n,T) ~ 3 • 10 9 cm n" 1 T 7 2 



(4) 



The strong Independence tells us that knowing the 
temperature is crucial. Note that n e and T e that enter 
here are those at the center of a Sweet-Parker reconnec- 
tion layer and are not known a priori. Therefore, we 
would like to cast condition (j4|) in an alternative form 
that would involve only the ambient plasma parameters, 
such as the far-upstream values of the plasma density, 
temperature, and magnetic field. Now, the cross-layer 
pressure balance (2n e T e = Bq/8tt), valid in the zero- 
guide-field case, provides us with one relationship be- 
tween the central n e and T e , and so, by itself, it is not 
sufficient. Indeed, this condition only tells us that the 
thermal pressure at the center of the layer needs to be 
raised to a certain level to balance the outside magnetic 
pressure, but it doesn't tell us whether this is achieved 
by increasing the density or the temperature. In order to 
break this degeneracy, we need to consider also the equa- 
tion of energy conservation. The logic of our model dic- 
tates that this analysis be done in the collisional Sweet- 
Parker regime. At the minimum, this analysis should 
include ohmic heating and heat advection and an esti- 
mate of various possible energy-loss mechanisms, such as 
radiation and electron thermal conduction. In particu- 
lar, it can be shown [28| that: (i) on the time that a 
fluid element spends inside the layer (the Alfven transit 
time t a = L/Va) — ohmic heating converts to heat just 
enough magnetic energy to raise T e to the level required 
by the pressure balance, without the need to increase 



the density substantially; (ii) for the solar coronal condi- 
tions, radiative losses are negligible on the Alfven time- 
scale; and (iii) the energy losses due to parallel electron 
thermal conduction can be neglected provided that the 
layer is collisional, in the sense of Eq. ((3|) ; (iv) the energy 
losses due to the perpendicular electron thermal conduc- 
tion are only marginally important at best. All these 
findings lead us to the conclusion that, once the condi- 
tion L > 40 A m f p is satisfied and hence the system is in 
the collisional Sweet-Parker regime, the energy equation 
can be regarded basically as a balance between ohmic 
heating and advection. As a result, the plasma density 
in the center of the layer should remain roughly compara- 
ble to the ambient coronal density, whereas the tempera- 
ture should increase dramatically. In fact, ohmic heating 
is enough to raise the temperature up to the "equipar- 
tition" level that depends only on the ambient density 
and upstream magnetic field Bq, and is insensitive to the 
ambient coronal temperature: 
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where -B1.5 = £?o/(30G). This estimate applies only as 
long as the resulting value is much higher than the am- 
bient temperature (typically of the order of 2 • 10 6 K for 
the solar corona), which means implies that the ambi- 
ent plasma-/? with respect to the reconnecting magnetic 
field should be <■ 1. Using this estimate, the collisionless 
reconnection condition can finally be written as [28L [2{| 



L < L c {n,B ) 



6- 10 9 cm n^ 3 Bf 5 



(6) 



In this form, the fast-reconnection condition is particu- 
larly useful because it involves only the global quantities 
characterizing a given reconnecting system: its global 
length L, the reconnecting component of the magnetic 
field So, and the ambient plasma density n e . 

Next, let us consider the strong-guide field case, B z ^> 
Bq, which is in fact more relevant to the problem of so- 
lar coronal heating. Although some of the arguments 
and results presented above have to be modified, concep- 
tually, they remain similar. In particular, the relevant 
collisionless-reconnection scale becomes the ion-acoustic 
gyro-radius, p s , calculated with the total magnetic field 
£?tot — B z (e.g., (13|). Correspondingly, the collisionless 
reconnection condition becomes 1311: 



Ssp < Ps 



a i Pe g 
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where, again, (3 e is based on the central n e and T e and 
on the upstream reconnecting field component Bq. Once 
again, all the quantities entering Eq. are to be es- 
timated in the collisional Sweet-Parker regime. To do 
this, first note that, in the strong guide field case one 
can no longer use the cross-layer pressure balance to 
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deduce /3 e ~ 1; this is because a relatively slight com- 
pression of the guide field can always ensure the pres- 
sure balance. Moreover, a strong guide field effectively 
makes the plasma incompressible, so that the central n e 
is equal to the ambient value. But this still leaves us 
with the task of evaluating the central electron temper- 
ature that is needed to determine the Spitzer resistivity. 
It turns out, however, that the above energy-balance ar- 
guments for a collisional Sweet-Parker layer still apply, 
at least qualitatively (28|. Therefore, the "equipartition" 
estimate for the central temperature, given by Eq. |(SJ), 
should still approximately hold; in particular, one should 
still have /3 e ~ 1. Then, one can repeat the procedure 
outlined above and derive the following approximate con- 
dition for the transition to fast collisionlcss rcconnection 
in the strong guide case [jjj ]: 



L < L c (n e ,B ,B z ) 



/mi f B 



6.10 9 cmn ro 3 ^. 5 (|l) .(8) 



Thus, the main effect of a strong guide field is to reduce 
the critical global length L c by a factor (B z /Bq) 2 3> 1; 
that is, the collisionlcss reconnection condition becomes 
harder to satisfy. Also, it is interesting to note that, for 
fixed values of n e and B Zl L c becomes very sensitive to 
the reconnecting field component: L c ~ Bq. 

Let us now discuss the implications of these findings for 
the solar corona. As long as flux emergence and braiding 
of coronal loops by photospheric footpoint motions con- 
tinue to generate current sheets in the corona, magnetic 
dissipation in these current sheets results in intermittent 
heating 0, Typical dimensions and field strengths 
of these current sheets are basically determined by the 
emerging magnetic structures and by the footpoint mo- 
tions. The main focus of this paper, then, is on what sets 
the typical level of the coronal plasma density and how it 
relates to the intermittent nature of energy release in the 
corona. My main point is that coronal heating should be 
viewed as a self-regulating process that keeps the corona 
marginally collisionless in the sense of Eqs. ((T|) — dHJ) (see 



HQ). 



As a first example of how this works, consider a coro- 
nal current sheet with some given fixed L, Bq, and B z . 
Resolving ([8]) with respect to n e , we can define a critical 
density, n c , below which reconnection switches from the 
slow collisional to the fast collisionless regime: 



n c {B z > B ) ~ 2 ■ 10 iu cm 



-3 R 4/3 .-1/3 
#1.5 L 9 
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Notice that this value is comparable to the typical den- 
sities observed in the active solar corona. I argue that 
this is just not a coincidence. Indeed, if at some initial 
time, the ambient density n e is higher than n c (L, Bq, B z ), 



then the current layer is collisional and reconnection is in 
the slow mode. Energy dissipation then is weak; hence, 
the plasma gradually cools through radiation and pre- 
cipitates to the surface. The density around the given 
current sheet drops and, at some point, becomes lower 
than n c . Then, the system switches to the fast collision- 
less regime and the rate of magnetic dissipation jumps. 
Next, there is an important positive feedback between 
coronal energy release and the density. A large part of 
the released energy is conducted to the surface, where 
it is deposited in a dense plasma. This leads to chro- 
mospheric evaporation along the post-reconnected mag- 
netic loops, filling them with a dense and hot plasma. 
The density rises and may now exceed n c . This will 
shut off any further reconnection (and hence heating) in- 
volving these loops until they again cool down, which 
occurs on a longer, radiative time-scale. Thus we see 
that, although highly intermittent and inhomogeneous, 
the corona is working to keep itself roughly at the critical 
density given by Eq. In this sense, coronal heating 
is self-regulating [291 ]. 

As a second example, consider a situation in which the 
initial density is relatively low, so that radiative cool- 
ing rate is much slower than the footpoint twisting rate. 
Then one can regard the density as constant between 
reconnection events and focus instead on the slow evo- 
lution of the reconnecting magnetic field, caused by the 
motion of the footpoints (similar to Refs. 0, [33J). Let us 
consider, for example, a flux tube, anchored on the solar 
surface at both ends, with a fixed strong axial (guide) 
field B z and a fixed volume V. The tube is comprised of 
smaller flux fibrils that are being wrapped around each 
other by the photospheric motions. Over time, this wrap- 
ping leads to the formation and strengthening of current 
layers. For simplicity, let us represent this process by a 
linear growth of the reconnecting field component of a 
single current sheet of a fixed length L: Bo(t) = ^tB z . 
Here, 7 parametrizes the rate of twisting. Let us now 
try to follow the evolution of this system. At first, Bq 
increases steadily in time, while the density stays con- 
stant. This continues until Bq reaches a critical value 
that depends on n e according to collisionlcss reconnec- 
tion condition (for fixed L and B z ): 



B c (n ei L,B z 



An (1 Ll /6 ntf Bl£ . (10) 



where B z _2 = £? Z /(100G). We assume here that Bo 
always stays well below B z . As soon as this critical 
value is reached, the system switches to the fast recon- 
nection regime and magnetic energy Bq/8it is rapidly 
dissipated. Importantly, part of this energy is not ra- 
diated promptly but is transported by parallel thermal 
conduction to the solar surface. This causes an evapora- 
tion episode adding new plasma to the flux tube under 
consideration. The amount of plasma added is roughly 
proportional to the energy released in a given event, 
Bq/8it ~ B^{n e )/8ir, which, according to (|10[) . is in turn 
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proportional to the density in the tube just before recon- 



nection: Sr 



Bi(n e ) 
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Now let us see what happens on a still longer (several 
hours) time-scale. As a consequence of the first fast re- 
connection event, the current sheet is promptly destroyed 
and Bq drops back to nearly zero. The field-line twist- 
ing, however, still continues, and so the process described 
above repeats. This time, however, the plasma density 
in the tube is higher, and hence the critical magnetic 
field B c is larger and takes longer time to reach. In par- 
ticular, taking the twisting rate 7 to be constant, the 
time between subsequent reconnection events scales as 
St = "f~ 1 B c (n e )/B z n X J 2 . Therefore, as long as the 
relative increase in density at each step is small, the long- 
term (t 3> St) evolution can be effectively described by 
the differential equation: 



dn e ^ Sn e (n e ) 
dt St(n e ) 



(11) 



and so n e (t) ~ t 2 . Correspondingly, the emission mea- 
sure of the tube increases as i 4 . 

This growth will continue until one of the following 
two effects intervene. First, as the density builds up, 
the critical value of Bq may become so large (a sizable 
fraction of B z ), that the equilibrium shape of the entire 
loop will be affected. The loop may then undergo the 
kink instability and become sigmoidal, which, with fur- 
ther twisting, may result in a large-scale eruption with a 
catastrophic energy release (a large flare). 

Alternatively, it may happen that the density just 
builds up gradually to a level large enough for radia- 
tive cooling between two subsequent reconnection events 
to become important. Indeed, as the density increases, 

the radiative emission measure increases as n 2 and the 

1/2 

time St between reconnection events as n e (see above). 
Ignoring for simplicity coronal temperature variations, 
the amount of thermal energy lost between two reconnec- 
tion events scales as n b J 2 , whereas the amount of ther- 
mal energy gained after each reconnection event is just 



proportional to B* 



At some point, the two 



will inevitably become comparable. Correspondingly, the 
amount of plasma drained due to the gradual radiative 
cooling will become equal to that pumped back up into 
the corona by each chromosphcric evaporation episode. 
Then, on some long time-scale (but still only as long as 7, 
L, and B z remain constant), the evolution of the system 
can be represented by repeated cycles that include fast 
reconnection events, followed by chromospheric evapo- 
ration episodes, followed by relatively long (~ 1 hour) 
periods during which the free magnetic energy builds up 
and the plasma gradually cools down. 
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